The Evolving Luminosity Function of Red Galaxies 
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ABSTRACT 

We trace the assembly history of red galaxies since z = 1, by measuring their evolving space 
density with the i3-band luminosity function. Our sample of 39599 red galaxies, selected from 
6.96 deg^ of imaging from the NOAO Deep Wide-Field and Spitzer IRAC Shallow surveys, is an 
order of magnitude larger, in size and volume, than comparable samples in the literature. We 
measure a higher space density of z ~ 0.9 red galaxies than some of the recent literature, in part 
because we account for the faint yet significant galaxy flux which falls outside of our photometric 
aperture. The -B-band luminosity density of red galaxies, which effectively measures the evolution 
of ^ L* galaxies, increases by only 36± 13% from z = to z = 1. If red galaxy stellar populations 
have faded by ~ 1.24 i?-band magnitudes since z — 1, the stellar mass contained within the red 
galaxy population has roughly doubled over the past 8 Gyr. This is consistent with star-forming 
galaxies being transformed into < L* red galaxies after a decline in their star formation rates. 
In contrast, the evolution of ~ 4L* red galaxies differs only slightly from a model with negligible 
z < 1 star formation and no galaxy mergers. If this model approximates the luminosity evolution 
of red galaxy stellar populations, then ~ 80% of the stellar mass contained within today's 4L* 
red galaxies was already in place at z = 0.7. While red galaxy mergers have been observed, such 
mergers do not produce rapid growth of AL* red galaxy stellar masses between z = 1 and the 
present day. 
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Red galaxies contain the majority of the stellar 
mass at low redshift (Hogg et al. 2002), so under- 
standing their formation and assembly is one of 
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the key goals of both observational and theoretical 
extragalactic astronomy. The stellar populations 
of these galaxies are dominated by an old com- 
ponent, with little ongoing star formation (e.g., 
Tinsley 1968). Observationally, it is unclear if this 
star formation occurred in situ at z > 1, or if 
these stars formed in lower mass galaxies which 
were later assembled via galaxy mergers. Sim- 
ulations of red galaxy evolution predict varying 
assembly and star formation histories, with some 
papers concluding that assembly takes place pri- 
marily at high redshift (e.g., Meza et al. 2003; 
Naab et al. 2005) while others suggest it contin- 
ues today (e.g., De Lucia et al. 2006). 

The color-magnitude diagram of low redshift 
galaxies provides several important clues about 
the assembly of red galaxies. At low redshift there 
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is a tight locus of red galaxies in color-magnitude 
space (e.g., Bower, Lucey, & Ellis 1992; Hogg et al. 
2004; Mcintosh et al. 2005), with the most lumi- 
nous red galaxies being slightly redder than less 
luminous red galaxies. This indicates that red 
galaxies of a given luminosity and metallicity have 
similar star formation histories, at least over the 
last few Gyr. A red galaxy locus is also observed 
within the 2; ~ 1 cluster and field galaxy popula- 
tions (van Dokkum et al. 2000; Blakeslee et al. 
2003; Bell et al. 2004; Willmer et al. 2005). The 
absence of very massive star-forming galaxies at 
low redshifts indicates that the most massive red 
galaxies are either formed at high redshift or they 
are assembled via mergers at z <1. As the colors 
of red galaxies are a function of luminosity, merg- 
ers of red galaxies should produce galaxies slightly 
blueward of the color-magnitude relation, and thus 
increase the scatter within the relation. 

Mergers of red galaxies, perhaps without signif- 
icant merger-triggered star formation (dry merg- 
ers), have been observed at low redshift (e.g., 
Lauer 1988; van Dokkum 2005), but it is unclear 
how large a role they play in galaxy formation. 
While there have been valiant attempts to mea- 
sure the galaxy merger rate with redshift (e.g., Le 
Fevre et al. 2000; Bell et al. 2006), there is debate 
about the selection function of merger candidates 
and the duration of observable galaxy mergers. At 
low redshift, measured rates of red galaxy stel- 
lar mass growth via mergers span from ~ 1% per 
Gyr (Masjedi et al. 2006) to ~ 10% per Gyr (van 
Dokkum 2005). 

The space density of red galaxies with redshift 
is a relatively dirc^ct measure of the galaxy assem- 
bly history, and lias fewer uncertainties than mea- 
surements of the galaxy merger rate. At the red- 
shifts where galaxy assembly takes place, the space 
density of galaxies as a function of stellar mass 
will evolve with redshift. While conceptually sim- 
ple, it is difficult to measure the galaxy space den- 
sity accurately, as z = 1 L* red galaxies are opti- 
cally faint and their strong spatial clustering (e.g.. 
Brown et al. 2003) results in significant cosmic 
variance. Despite these difficulties, several groups 
have used measurements of the space density of 
galaxies with redshift to constrain red galaxy evo- 
lution and assembly (e.g., Lilly et al. 1995; Lin 
et al. 1999; Chen et al. 2003; Bell et al. 2004; 
Bundy et al. 2005; Willmer et al. 2005; Wake 



et al. 2006). Faber et al. (2005) provides a useful 
summary of prior studies, and describes their vary- 
ing conclusions. The prior literature provides sev- 
eral plausible scenarios for z < I red galaxy evolu- 
tion, including passive evolution, assembly via dry 
mergers, formation from fading blue galaxies, or a 
combination of the above. 

This is the first in a series of papers which will 
discuss the assembly and evolution of z < 1 red 
galaxies detected by the multiwavelength imaging 
surveys of the entire Bootes field. Other papers 
in this series will measure the stellar mass func- 
tion, spatial clustering, and AGN content of red 
galaxies. In this paper we present the evolving 
B-band luminosity function of 0.2 < 2; < 1.0 red 
galaxies, and discuss the resulting implications for 
red galaxy evolution and assembly. While the 
rest-frame B-band is not ideal for tracing stel- 
lar mass, it has the advantages of being within 
the well studied optical wavelength range and re- 
maining within the observed BwRI passbands of 
the NOAO Deep Wide-Field Survey (NDWFS) for 
redshifts of z < 1. 

The structure of this paper is as follows. We 
provide an overview of the NDWFS and Spitzer 
IRAC Shallow Survey in §2. We discuss our 
catalogs, photometry, photometric redshifts, and 
galaxy rest-frame properties in §3. In §4 we de- 
scribe the selection of our red galaxy sample. In 
§5 we present rest-frame _B-band luminosity func- 
tions of red galaxies and compare our measure- 
ments with simple galaxy evolution models. In 
§6 we compare our results with recent red galaxy 
surveys, and discuss discrepancies between various 
studies. Our principal results and conclusions are 
summarized in §7. Throughout this paper we use 
Vega based magnitudes and adopt a cosmology of 
nm = 0.3, r^A = 0.7, Ho = 70 km s"^ Mpc"\ 
h = 0.7, and erg = 0.8. 

2. THE SURVEYS 

2.1. The NOAO Deep Wide-Field Survey 

The NOAO Deep Wide-Field Survey (NDWFS) 
is an optical {BwRI) and near-infrared {K) imag- 
ing survey of two « 9.3 deg^ fields with telescopes 
of the National Optical Astronomy Observatory 
(Jannuzi & Dey 1999). A thorough description of 
the observing strategy and data reduction will be 
provided by B. T. Jannuzi et al. and A. Dey et 
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al. (both in preparation). Wo utilize the third 
NDWFS data release^ of optical imaging with the 
MOSAIC-I camera on the Kitt Peak 4-m tele- 
scope. To obtain accurate optical colors with fixed 
aperture photometry across the Bootes field, we 
have smoothed copies of the released images so the 
stellar Point Spread Function (PSF) is a Moffat 
profile with a full width at half maximum of 1.35" 
and 13 = 2.5. The subfield NDWFSJ1428 + 3531 
(each subfield is roughly one MOSAIC-I field-of- 
view) has poor seeing in the Bw and i?-bands, and 
has been excluded from this study. 

2.2. The IRAC Shallow Survey 

The IRAC instrument (Fazio et al. 2004) on 
the Spitzer Space Telescope provides simultaneous 
broad-band images at 3.6, 4.5, 5.8, and 8.0 /./,m. 
The IRAC Shallow Survey, a guaranteed-time ob- 
servation program of the IRAC instrument team, 
covers 8.5 deg^ of Bootes with three or more 
30 second exposures per position. Eisenhardt 
et al. (2004) present an overview of the survey 
design, reduction, calibration, and preliminary re- 
sults. Despite the short integration times, the 
IRAC Shallow Survey easily detects 2; ~ 1.4 clus- 
ter galaxies (Stanford et al. 2005; Elston et al. 
2006; Brodwin et al. 2006). In this paper we uti- 
lize the 3.6 and 4.5 /zm imaging to remove contam- 
inants (e.g., stars, quasars) from our galaxy sam- 
ple and for empirical photometric redshifts. The 
observed / — [3.6] color of ^; < 1 red galaxies is 
a superb redshift indicator, due to the redshift- 
ing of the 1.6 /xm H~ stellar opacity feature (e.g., 
Simpson & Eisenhardt 1999; Sawicki 2002). 

3. THE OBJECT CATALOG 
3.1. OBJECT DETECTION 

We detected sources using SExtractor 2.3.2 
(Bertin & Arnouts 1996), run in single- image 
mode on the /-band images of the NDWFS third 
data release. The NDWFS Bootes field comprises 
of 27 optical subfields, each of which are 35' x 35' 
in size, that are are designed to have overlaps of 
several arcminutes. Individual objects can be de- 
tected in multiple subfields, which is particularly 

^Available from the NOAO Science Archive at 
http://www.aj:chive.noao.edu/ndwfs/ 



useful for photometric calibration, though the ef- 
fective exposure time generally decreases towards 
subfield edges. For this paper, we only include ob- 
jects which are detected within nominal subfield 
boundaries, which we define to have overlaps of 
only tens of arcseconds. In these small overlap re- 
gions, we search for objects with /-band detections 
in multiple subfields and retain the detection with 
the highest quality data (i.e., without bad pixels 
or with the longest exposure time). To minimize 
the number of faint spurious galaxies in our cata- 
log, we exclude regions surrounding very extended 
galaxies and saturated stars. We also exclude re- 
gions that do not have good coverage from both 
the NDWFS and the IRAC Shallow Survey. The 
final sample area is 6.96 deg^ over a 2.7° x 3.3° 
field-of-view. 

3.2. PHOTOMETRY 

We used our own code to measure aperture 
photometry for each object in the optical and 
IRAC passbands. To reduce contaminating flux 
from neighboring objects, we used SExtractor seg- 
mentation maps to exclude pixels associated with 
neighboring objects detected in any of the three 
optical bands. We corrected the aperture pho- 
tometry for missing pixels (from the segmentation 
maps or bad pixels) by using the mean flux per 
pixel measured in a series on 0.5" wide annuli sur- 
rounding each object position. Accurate la un- 
certainties for the photometry were determined by 
measuring photometry at ~ 100 positions within 
2' of the object position and finding the uncer- 
tainty which encompassed 68.7% of the measure- 
ments. We also use the median of the fluxes mea- 
sured at these positions to subtract small errors 
which could be present in our sky background es- 
timate. Typical random uncertainties for /-band 
4" diameter aperture photometry are 0.03 mag- 
nitudes at / = 21, increasing to 0.2 magnitudes 
at / = 23. Typical random uncertainties for 
3.6 fim 4" aperture photometry are 0.02 magni- 
tudes at [3.6] = 16, increasing to 0.2 magnitudes 
at [3.6] = 19. 

To verify the accuracy of our uncertainties and 
to search for systematic errors, we added artiflcial 
galaxies to our data, recovered them with SExtrac- 
tor, and measured their photometry. The artificial 
galaxies have de Vaucouleurs (1948) profiles trun- 
cated at seven half-light radii (though our results 
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Fig. 1. — The offset between measured and total photometry for 4" aperture photometry (left) and SExtrac- 
tor's MAG_AUTO (right). We only plot galaxies over the apparent magnitude range where the completeness 
is 50% or more, and we treat non-detections as measurements of / = 99. The predicted offsets due to flux 
outside the 4" aperture are shown with arrows, and these are consistent with our measurements. SExtractor's 
MAG_AUTO photometry is closer to total magnitudes than 4" aperture photometry, but has systematics as 
a function of both half-light radius and apparent magnitude. Luminosity functions using MAG_AUTO for 
total magnitudes may have systematics as a function of both luminosity and redshift. 



are not particularly sensitive to the truncation ra- 
dius), which were then convolved with our PSF. 
Figure 1 plots the offset between total and mea- 
sured luminosities for these galaxies. This offset 
is due to flux outside outside of our 4" diameter 
aperture, and is a predictable function of half-light 
radius. The observed half-light radius is a func- 
tion of both galaxy luminosity and redshift, and 
we discuss corrections from 4" aperture photome- 
try to total magnitudes in §3.4. After accounting 
for this predicted offset, we find that the differ- 
ence between the input and measured photometry 
is within the quoted la uncertainties ~ 70% of the 
time. 

It is useful to compare our photometry with 
SExtractor's MAG_AUTO, which is often used 
to measure "total" magnitudes in galaxy surveys 
(e.g., Bell et al. 2004; Zucca et al. 2005). While 
SExtractor is an extremely useful tool for source 
detection, photometry and classification, it can- 
not be expected to provide perfect measurements 
of all objects in all surveys. Using the same ar- 
tificial galaxies as described above, we find 80% 
of the MAG_AUTO measurements differ from the 



total magnitude by more than the quoted Icr un- 
certainty. This is due to SExtractor's assumption 
that the sky background is Gaussian random noise 
without source confusion, which is only a valid ap- 
proximation for some imaging data. 

We show in Figure 1 that MAG_AUTO is closer 
to a total magnitude than 4" aperture photome- 
try, though it has systematic errors which are a 
function of apparent magnitude. In addition, its 
random errors are 40% larger than 4" aperture 
photometry. The systematic error with magnitude 
may be due to the elliptical MAG_AUTO aperture 
being defined using the second-order image mo- 
ments of object pixels above an isophotal thresh- 
old. It is therefore plausible that the MAG_AUTO 
aperture is slightly too small for relatively ex- 
tended faint galaxies. Luminosity functions using 
MAG_AUTO for total magnitudes may have sys- 
tematics of 0.2 magnitudes at faint magnitudes, 
though the exact size of this offset will depend 
both on the imaging data and user defined SEx- 
tractor parameters. Accurate luminosity functions 
require corrections for galaxy flux which falls be- 
yond the photometric aperture. 
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3.3. PHOTOMETRIC REDSHIFTS 

Wc determined redsliifts for our galaxies us- 
ing the empirical ANNz photometric redshift code 
(Firth, Lahav, & Somerville 2003; Collister & La- 
hav 2004). When a sufficiently large training set 
is available, empirical photometric redshifts are 
as precise or better than other techniques (e.g., 
Csabai et al. 2003; Brodwin et al. 2006). This 
approach does restrict us to z < 1, where thou- 
sands of spectroscopic redshifts are available in 
Bootes to calibrate the photometric redshifts. 
ANNz uses artificial neural networks to determine 
the relationship between measured galaxy proper- 
ties and redshift. It does not use any prior assump- 
tions about the shape of galaxy spectral energy 
distributions (SEDs), though it does assume the 
relationship between galaxy properties and red- 
shift is a relatively smooth function. ANNz works 
best when the training set is a large representative 
subset of the science galaxy sample. 

The basis of our training set is galaxies with 
spectroscopic redshifts in the Bootes field. The on- 
going AGN and Galaxy Evolution Survey (AGES, 
C. S. Kochanek et al. in preparation) of Bootes has 
obtained spectroscopic redshifts of ~ 16000 / < 20 
galaxies, using the MMT's Hectospec multiobject 
spectrograph (Fabricant et al. 2005). At fainter 
magnitudes, most of the spectroscopic redshifts 
are from various spectroscopic campaigns with the 
W. M. Keck, Gemini, and Kitt Peak observatories, 
which have jointly resulted in redshifts for several 
hundred galaxies. 

We trained and measured photometric redshifts 
using the 4" aperture photometry and the 2nd 
order moments of the /-band light distribution. 
For training and determining photometric red- 
shifts we used asinh magnitudes (Lupton, Gunn, 
& Szalay 1999) rather than fluxes or conventional 
2.51og]^Q magnitudes. ANNz does not produce ac- 
curate photometric redshifts with fluxes that span 
a vast dynamic range, while conventional magni- 
tudes only measure positive fluxes. We only use 
asinh magnitudes when estimating ANNz photo- 
metric redshifts and use conventional magnitudes 
throughout the remainder of the paper. The sec- 
ond order image moments are usehil for reduc- 
ing the number of bright galaxy photometric red- 
shifts with large errors. In particular, they re- 
duce photometric redshift outliers produced by 



low redshift edge-on spirals whose dust lanes pro- 
duce colors similar to higher redshift objects. As 
we do not want the shape measurements to be 
a function of /-band depth, which varies slightly 
across Bootes, we measure the 2nd order moments 
with pixels above an /-band surface brightness of 
22.5 mag arcsec"^. 

We applied several restrictions to the training 
set to improve the photometric redshifts of 2 < 1 
red galaxies. We restricted the training set to 
z < 1.5 spectroscopic galaxies, so ^; > 1.5 galaxies 
did not perturb the ANNz photometric redshift 
solution for z < 1.5 galaxies. While we include 
red galaxies with X-ray and radio counterparts in 
our final sample, these objects are overrepresented 
in our spectroscopic samples and a small fraction 
of these sources could have unusual colors due to 
the contribution of the AGN. We have therefore 
excluded objects with counterparts in the Chan- 
dra XBootes survey (Kenter et al. 2005; Murray 
et al. 2005; Brand et al. 2006) or the FIRST ra- 
dio survey (Becker, White, & Helfand 1995) from 
the photometric redshift training set. The small 
number of galaxies with unusual apparent colors 
which remained in the training set were removed 
with apparent color cuts which are a function of 
spectroscopic redshift. 

At z > 0.6, there are only 280 red and blue 
galaxies fainter than / = 20.5 with spectroscopic 
redshifts in the training set, which is less than ideal 
for training ANNz. As the shape of the galaxy 
color-locus at any given redshift is relatively sim- 
ple and smooth, we were able to increase the size 
of the training set by adding interpolated objects 
to the training set. We did this when we found 
galaxies within 0.02 in redshift and within 0.5 in 
B\Y — R color of each other. As red galaxies at 
z < 1 follow relatively tight color-magnitude (e.g., 
Bower et al. 1992; Bell et al. 2004; Mcintosh et al. 
2005) and size-luminosity relations (e.g., Bender, 
Burstein, & Faber 1992; Shen et al. 2003), we 
were able to extrapolate the training set to fainter 
magnitudes by making faint copies of bright ob- 
jects with slightly altered photometry and image 
size parameters. As the colors and sizes of galaxies 
are a very strong function of redshift, approxima- 
tions of the color-magnitude and size-luminosity 
relations at a given redshift are sufficient for the 
photometric redshift training set. Figure 2 plots 
the color-color and color-magnitude diagrams of 



5 





Fig. 2. — Apparent color-color and color-magnitude diagrams for photometric redshift training set galaxies 
with spectroscopic redshifts 0.78 < z < 0.82. Bold symbols denote real galaxies while artificial objects 
produced by interpolating in color or extrapolating in magnitude arc shown with dots. As the observed 
colors of galaxies are a strong function of redshift and a weak function of luminosity, approximations of the 
color-color and color-magnitude relations can be used to extrapolate the training set to faint magnitudes. 



training set galaxies at z = 0.80, along with inter- 
polated and extrapolated objects. 

In Figure 3 we plot the photometric and spec- 
troscopic redshifts of 4314 objects which meet our 
red galaxy selection criteria (which we discuss in 
§4), including radio and X-ray sources. We also 
provide estimates of photometric redshift uncer- 
tainties for a series of redshift and apparent mag- 
nitude bins in Table 1. The uncertainties were 
determined using the relevant percentiles, and do 
not assume a Gaussian distribution of photometric 
redshift errors. As the observed range of galaxy 
colors broadens with increasing apparent magni- 
tude, at the very least due to photometric errors, 
the random errors of the photometric redshifts also 
increase with apparent magnitude. Measured sys- 
tematic errors are much smaller than the random 
errors listed in Table 1. The accuracy of our red 
galaxy photometric redshifts are comparable to 
the best broad-band photometric redshifts avail- 
able in the literature (e.g., Mobasher et al. 2004). 
The la uncertainties of our photometric; redshifts 
are ~ 0.1 in redshift at / = 22, and this decreases 
to ~ 0.03 at 7 = 19.5. 
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Fig. 3. — Comparison of spectroscopic and photo- 
metric redshifts for the red galaxy sample. Sym- 
bol greyscale is a function of apparent magnitude, 
with the small number of 7 > 22 galaxies with 
spectroscopic redshifts having large black symbols. 
A total of 4314 objects are plotted, and only 0.3% 
of these objects have errors of > 0.2 in redshift. 
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3.4. REST-FRAME PROPERTIES 

To measure the rest-frame absolute magnitudes 
and colors of NDWFS galaxies, we used maxi- 
mum likelihood fits of Bruzual & Chariot (2003) 
SED models to the BwRI 4" diameter aperture 
photometry. Throughout this paper we use so- 
lar metallicity models with a Salpeter (1955) ini- 
tial mass function, a formation redshift of z = 4, 
and exponentially declining star formation rates. 
These models provide a reasonable approximation 
of the observed SEDs of red galaxies. Unlike em- 
pirical templates, the Bruzual & Chariot (2003) 
models allow us to estimate stellar masses and star 
formation rates, which we will use in future pa- 
pers (K. Brand et al. and A. Dey et al., both in 
preparation). For this paper, we will only use the 
rest-frame optical colors and absolute magnitudes. 

Other stellar population models can be fitted 
to red galaxies, though changing models has little 
impact on our B-band hmiinosity fmictions and 
our conclusions. Shifting the formation redshift 
results in different r models being fitted to the red 
galaxies, but the model i?-band luminosity evolu- 
tion of red galaxy stellar populations at z < 1 
changes by 0.1 magnitudes or less. Changing from 
a Salpeter (1955) to a Chabrier (2003) initial mass 
function decreases galaxy stellar masses by 25%, 
but has little impact on the predicted evolution 
of galaxy optical colors and Imiiinosities. Models 
with low metallicities are bluer than the reddest 
galaxies, while models with 1.5 times Solar metal- 
licity are offset from observed BwRI galaxy loci 
from z = 0.2 to z 1.0. 

As with all current galaxy SED templates and 
models, there are small systematic errors. The r 
models do not perfectly match the evolving rest- 
frame U —V colors of red galaxies at all redshifts, 
and the t models overestimate the apparent R — I 
colors of 0.4 ^ z < 0.7 red galaxies by ~ 0.05 
magnitudes. As discussed in §6, we expect these 
systematics to have a modest impact on our lumi- 
nosity functions. 

The 4" aperture photometry captures 86% or 
less of the total flux. We corrected for the flux 

outside this aperture by assuming galaxies within 
our sample have truncated de Vaucouleurs (1948) 



profiles, 



/ a exp 
— exp 



-7.6695 I — 



1/4- 



-7.6695 X 7^/^ 



(1) 



at r < 7re, where the half-light radius Ve is a func- 
tion of S-band absolute magnitude and redshift. 
At 2; = 0, we use the size- luminosity relation of 
SDSS early-type galaxies (Shen et al. 2003) and 
assume Bvega — ^ab = 1.32 (Fukugita, Shimasaku, 
& Ichikawa 1995), so 

log[re(h-ikpc)] = 1.0 -0.26(MB-51ogh-F 21.81). 

(2) 

The la dispersion around this relation is only 30% 
for the most luminous red galaxies (Shen et al. 
2003), while the offsets between 4" and total mag- 
nitudes in Figure 1 are a relatively weak function 
of half-light radius. This simple relation should 
therefore provide accurate corrections between 4" 
aperture and total photometry. 

The size-luminosity relation must change with 
redshift due to the evolution of red galaxy stel- 
lar populations. We assume the size- luminosity 
relation undergoes pure luminosity evolution de- 
scribed by a Bruzual & Chariot (2003) r = 0.6 Gyr 
stellar synthesis model with a formation redshift 
of 2; = 4. This model fades by 1.24 B-band mag- 
nitudes between z = 1 and z = 0. As shown in 
Figure 4, the r = 0.6 Gyr model approximates 
the color evolution of My — 5 log h < —21 — z 
red galaxies. The luminosity evolution of the fun- 
damental plane for Mb — 5 log h ~ —20 — 1.2z 
early- type galaxies (e.g., van Dokkum & Stanford 
2003; Treu et al. 2005) is also comparable to this 
model, though less luminous galaxies probably ex- 
hibit more rapid luminosity function (Treu et al. 
2005; van der Wei et al. 2005). Stellar population 
synthesis models with r < 1 Gyr have less than 2% 
of their total star formation at z < 1. While this 
is a tiny amount of star formation, it does result 
in more rest-frame color and luminosity evolution 
than passive evolution models without any z < 1 
star formation. Changing the value of t by 0.2 Gyr 
alters the model z = 1 rest- frame U — V colors and 
B-band luminosity evolution by ~ 0.1 magnitudes. 
After accounting for the evolving size-luminosity 
relation and our 1.35" point spread function, we 
find the 4" diameter aperture captures 75% of the 
fiux of a Mb — 5 log h = —21 galaxy at z = 0.9. 
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Fig. 4. — The evolving rest-frame U— V colors of red and blue galaxies brighter than My — 5 log h = —21 — 2;. 
For comparison, the U — V colors of My — 5 log h = —21 red galaxies in Virgo (Bower et al. 1992) and 
z = 0.05 clusters (Mcintosh et al. 2005) are also shown. The locus of galaxies which lie along the red galaxy 
color-magnitude relation is evident. The evolution of the color- magnitude relation is well approximated by 
a Bruzual & Chariot (2003) r = 0.6 Gyr model (solid red line), which has an exponentially declining star 
formation rate, no mergers, and little ongoing star formation at z < 1. 
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Fig. 5. — The observed and predicted offset between 4" and 6" aperture photometry for red galaxies. For the 
vast majority of z > 0.4 red galaxies the measured offset is within 0.15 magnitudes of the value expected for 
point sources (0.11 magnitudes). The median of the observed offsets is shown by a purple dashed line in both 
panels. The SDSS size-luminosity relation (Shen et al. 2003) without luminosity evolution (red dash-dot 
line) overestimates the sizes of red galaxies. The SDSS size-luminosity relation with luminosity evolution 
described by a r = 0.6 Gyr Bruzual & Chariot (2003) model (blue solid line) approximates the data well. For 
a fixed apparent /-band magnitude, the observed half-light radius of red galaxies increases with increasing 
redshift. Estimators of total galaxy magnitudes may therefore exhibit larger systematic errors at z = 0.9 
than z ~ 0.5. 
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Red galaxies do not exclusively have de Vau- 
couleurs (1948) profiles (e.g., Blanton et al. 2003; 
Driver et al. 2006), and this could affect our esti- 
mates of galaxy total magnitudes. To investigate 
this, we determined the offsets between 4" aper- 
ture and total magnitudes for Sersic (1968) profiles 
with indices between 2 and 6. For galaxies with 
half-light radii of 1", the offset differs from that de- 
termined with a de Vaucouleurs (1948) profile by 
0.03 magnitudes. Changing the truncation radius 
from 7 to 10 half-light radii only changes the off- 
set between 4" aperture and total magnitudes by 
^0.06 magnitudes for galaxies with ^ 1". This 
offset can be increased by using larger truncation 
radii, though this results in a systematic difference 
between the model and measured offset of 4" and 
6" aperture photometry for red galaxies. As the 
expected evolution of red galaxy stellar popula- 
tions is on the order of 1.24 B-band magnitudes 
per unit redshift, small errors in our corrections 
from 4" aperture to total magnitudes should have 
little impact upon our results and conclusions. 

To verify the accuracy of our model, we com- 
pared the measured and predicted difference be- 
tween 4" and 6" aperture photometry for red 
galaxies. As shown in Figure 5, our simple model 
provides a good approximation of what is observed 
in our red galaxy sample. Figure 5 also illustrates 
the absolute magnitude dependence of our model. 
This dependence results in the aperture correc- 
tions at fixed apparent /-band magnitude increas- 
ing with redshift. It is therefore plausible that 
an estimator of total magnitudes could work well 
for the majority of / = 21 galaxies, which are at 
z < 0.7, but have systematic errors for / = 21 
galaxies at z = 0.9. Accounting for flux out- 
side the aperture is crucial for measuring the total 
magnitudes of 2; < 1 galaxies. 

4. THE RED GALAXY SAMPLE 

The rest-frame distribution of galaxy colors is 
bimodal (e.g., Hogg et al. 2004), and selection cri- 
teria for red galaxies typically fall near the mini- 
mum between the red and blue galaxy populations 
(e.g., Madgwick et al. 2002; Bell et al. 2004; Faber 
et al. 2005). We apply a similar approach for this 
work, and use the following rest-frame color selec- 
tion criterion; 



U-V > 1.40-0.25 

-0.08 X (My - 5 log /i -h 20.0) 

-0.42 X (z-0.05) 

-H0.07 X (z - 0.05)2. (^3-) 

Our criterion selects galaxies with rest-frame U — 
V colors within 0.25 magnitudes of the evolving 
color-magnitude relation of red galaxies. This cri- 
terion allows comparison with the recent literature 
and is very similar, though not identical, to the 
criterion of Bell et al. (2004). 

Our selection criterion is plotted in Figure 6 
along with the observed color distribution of red 
galaxies. Our selection criterion has slightly more 
tilt than the observed color-magnitude relation 
of 0.2 < z < 0.4 red galaxies. If we use 8" 
aperture photometry instead of 4" aperture pho- 
tometry, the observed color-magnitude relation of 
0.2 < z < 0.4 red galaxies moves blueward and 
steepens. However, as this has little impact upon 
the measured luminosity function, we continue to 
use 4" aperture photometry, which has smaller 
random uncertainties than 8" aperture photom- 
etry. 

The measured space density of the most lumi- 
nous red galaxies does not depend on the details 
of our selection criterion, as these galaxies mostly 
lie along the color-magnitude relation. We can 
therefore easily compare our luminosity functions 
for Mb — 5 log h < —21 red galaxies with those of 
2dFGRS, SDSS, COMBO-17, and DEEP2 (Madg- 
wick et al. 2002; Bell et al. 2004; Faber et al. 
2005; Willmer et al. 2005; Blanton 2005), all of 
which use slightly different galaxy selection crite- 
ria. The fraction of blue galaxies increases with 
decreasing luminosity, so the measured space den- 
sity of < L* red galaxies is a stronger function of 
the red galaxy selection criteria. If we shift our 
criterion blueward by 0.10 magnitudes our mea- 
sured space density of Mb — 5 log h ^ —20 red 
galaxies increases by ~ 25%. 

We limit the absolute magnitude range in each 
of our redshift bins so we can determine accurate 
redshifts and have a highly complete sample. Al- 
most all of the galaxies in our flnal sample are de- 
tected in Bw, R, I and the two IRAC bands. To 
evaluate the completeness of our sample, we added 
artificial galaxies to copies of the /-band images 
and recovered them with SExtractor. The arti- 
ficial galaxies had de Vaucouleurs (1948) profiles 
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Fig. 6. — Color-magnitude diagrams for red galax- 
ies. Our selection criterion is shown with the 
dashed blue line and is 0.25 magnitudes blue- 
ward of the model color-magnitude relation (solid 
red line). The color- magnitude relation moves to- 
wards bluer rest-frame colors with increasing red- 
shift. There is a difference between the model 
and observed color-magnitude relation slope at 
0.2 < z < 0.4, due to our use of 4" aperture pho- 
tometry, which measures the cores of low redshift 
galaxies. As discussed in §3.4, this has negligible 
impact on our luminosity functions and conclu- 
sions. 

truncated at 7 effective radii, and we measure the 
completeness for a range of effective radii and ap- 
parent magnitudes. As illustrated in Figure 7, our 
catalogs are more than 85% complete for / < 23.5 
galaxies with half-light radii of 0.5" or less. Our 
measurements of the luminosity function, which 
we discuss in §5, include small corrections for this 
incompleteness. 

4.1. ADDITIONAL SELECTION CRI- 
TERIA 

Contaminants cannot be rejected from photo- 
metric redshift surveys on the basis of their spec- 
tra, so these surveys can be more susceptible 
to contaminants than comparable spectroscopic 
surveys. This contamination can include stars. 



Fig. 7. — The completeness of the red galaxy sam- 
ple as a function of /-band apparent magnitude. 
We measured the completeness by adding artificial 
galaxies to copies of the data and recovering them 
with SExtractor. As our incompleteness is dom- 
inated by confusion with brighter (and generally 
lower redshift) sources rather than fluctuations in 
the sky background, our completeness is a weak 
function of half-light radius. The completeness is 
higher than 85% over our sample magnitude range, 
and is well approximated by 1 — 0.05(7 — 21) for 
21.0 < / < 23.5 galaxies. 

quasars, and galaxies with large photometric red- 
shift errors. We therefore apply apparent color 
and morphology cuts to minimize contamination 
while retaining a highly complete sample of red 
galaxies. 

We have applied apparent color cuts to exclude 
objects whose colors differ from those expected for 
0.2 < z < 1.0 red galaxies. These color cuts are 
designed to remove most stars, quasars and z ^ 1 
galaxies from our sample. Our cuts may exclude 
a small percentage of red galaxies with unusu- 
ally large photometric errors, but as these galax- 
ies could also have large redshift and luminosity 
errors, it is preferable to exclude them from our 
measurement of the luminosity function. As the 
distribution of red galaxy apparent colors does not 
form a simple shape in color space, we use a total 
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R-l l-[3.6] 



Fig. 8. — The apparent color distributions of Mb — 5\ogh < —19.5 red galaxies in Bootes. Red lines 
denote apparent color cuts, which we use to exclude stars, quasars, blue galaxies, z ^ 1 galaxies, and gross 
photometric outliers. For faint objects with highly uncertain apparent colors (e.g., Bw > 26.5 galaxies), we 
use upper or lower limits for the color when applying the color cuts. The locus of 18 < / < 19 stars and 
quasars is shown in blue in the right-hand panel, and these objects are largely excluded by our color cuts. 
The red galaxy color-magnitude relation produces an over-density in color space that is clearly evident. 



17 color cuts to exclude objects from the sample. 
All but two of these cuts are plotted in Figure 8, 
while a full list is provided in Table 2. 

Our color cuts remove 10767 objects from the 
sample, reducing it to 39866 objects. As we apply 
the color cuts before applying the morphology cri- 
terion, most of the objects rejected from the sam- 
ple have the colors and morphologies of stars and 
quasars. Unlike 0.2 < z < 1.0 red galaxies, main 
sequence stars have apparent colors which satisfy 
R~ I > 0.65 X (/ - [3.6] - 1) or / - [3.6] < 2.2, 
and 8775 of the objects excluded by our color cuts 
satisfy these two criteria. Another 749 of the ex- 
cluded objects satisfy [3.6] — [4.5] > 0.6, and are 
probably quasars. The measured photometry of 
the remaining 1223 objects excluded by our color 
cuts differs somewhat from main sequence stars, 
most quasars, and 0.2 < z < 1.0 red galaxies. 
Spectroscopic redshifts are available for 23 of these 
objects, and 13 are 0.2 < z < 1.0 galaxies while 
the remainder are stars and quasars. It is therefore 



plausible that ~ 700 red galaxies are excluded by 
our color cuts, though this is only 2% of our final 
sample. 

To further reduce contamination by compact 
objects, we also apply a morphology criterion. 
We exclude objects if the difference between their 
2" and 4" /-band aperture photometry is 0.20 
magnitudes below the expectation from the size- 
luminosity relation discussed in §3.4. We expect 
many of the 287 objects rejected by this criterion 
to be blends of galaxies with stars and quasars. 
The morphology and apparent color cuts could 
plausibly exclude ~ 10'^ red galaxies from our fi- 
nal sample. This reduces sample size by ~ 2.5%, 
but as our uncertainties from cosmic variance are 
on the order of 10% (§5), the impact upon our 
results and conclusions is insignificant. The fi- 
nal 0.2 < z < 1.0 red galaxy sample contains 
39599 galaxies brighter than / — 23.5, and number 
counts as a function of both redshift and apparent 
magnitude are summarized in Table 3. 
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5. THE RED GALAXY LUMINOSITY 
FUNCTION 

We measured the red galaxy luminosity func- 
tion in four redshift slices between z — 0.2 and 
z — 1.0. We used both the non-parametric 1/Vmax 
technique (Schmidt 1968) and maximum likeli- 
hood fits (e.g., Marshall et al. 1983) of Schechter 
(1976) functions; 

^{M)dM = O.41nlOx0* (^-^^ exp (^^^ dM, 

(4) 

where L is the galaxy luminosity while (f>* , L* , and 
a are constants. Like most luminosity function 
papers, we use M* rather than L*, where M — 
M* = -2.51og(L/L*). 

The uncertainties of the luminosity function 
are dominated by large-scale structure rather than 
Poisson counting statistics. As galaxies with dif- 
ferent luminosities can occur within the same 
large-scale structures, the data-points in binned 
luminosity functions are not independent of each 
other. We have evaluated the uncertainties of the 
luminosity function using both subsamples of the 
Bootes field and the galaxy angular correlation 
function. 

Subsamples are conceptually simple but under- 
estimate the uncertainties, as an individual large- 
scale structure may span 2 subsamples of the data. 
For this reason, we only use thirteen 0.5 deg^ 
subsamples rather than many smaller subsamples. 
For our Schechter function fits we evaluate the 
luminosity function for each subsample using the 
method of Marshall et al. (1983) and use the stan- 
dard deviation of the fitted parameters (e.g., M*) 
divided by -y/TS to estimate uncertainties. Lumi- 
nosity functions for the whole Bootes field and the 
thirteen subsamples are shown in Figure 9. The 
subample luminosity functions can differ from the 
luminosity function for the entire field by as much 
as 50%. While individual galaxy clusters are ev- 
ident in Figure 10, these contain but a fraction 
of all red galaxies and the variations between dif- 
ferent subsamples are almost certainly caused by 
galaxies residing within larger structures. 

As the angular correlation function of galaxies 
does not equal zero on scales of ~ 1°, we expect 
subsamples to underestimate the uncertainties for 
(j>* and the luminosity density, Jb- We do not cal- 




Fig. 9. — Red galaxy luminosity functions for the 
entire Bootes field and our thirteen subsamples. 
The solid lines are maximum likelihood Schechter 
function fits to the data while the symbols are 
^/Vmax estimates of the luminosity function. The 
luminosity functions for each 0.5 deg^ subsample 
are shown in a different greyscale while the lu- 
minosity functions for the entire Bootes field are 
shown in color. While individual galaxy clusters 
are evident in Figure 10, these contain only a frac- 
tion of all red galaxies and the variations between 
different subsamples are almost certainly caused 
by galaxies residing within larger structures. 

culate uncertainties for Al* and a using galaxy 
clustering, as this requires additional information 
including details of how the shape of the lumi- 
nosity function varies with galaxy density. The 
expected variance of the number counts in a field 
is given by 

(5) 

(Groth & Peebles 1977; Efstathiou et al. 1991) 
where i^{0) is the angular correlation function, 
is the angle separating solid angle elements dfli 
and dil2, and il is the area of the field. We as- 
sume i-u{0) is a power-law with index 1 — 7, and 
use power-law fits to the angular correlation func- 
tions from M. J. L Brown et al. (in preparation). 
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Fig. 10. — The sky distribution of red galaxies in the 2.7° x 3.3° Bootes field. In the electronic edition, 
0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8 and 0.8 < z < 1.0 objects are denoted with light blue, dark 
blue, green and red symbols respectively. Grey regions denote areas excluded from the sample, including 
saturated stars and areas without both NDWFS and Spitzer coverage. Individual structures are clearly 
evident, and 0.4 < z < 0.6 galaxy clusters are not uniformly distributed across the Bootes field. At z > 0.6, 
the distribution of galaxies is more uniform, though individual structures can still be seen. Of particular 
note is a - SOh'^Mpc long z ~ 0.93 structure at 14''35'"30% + 34°50'. 



For the Bootes field the sample variance is ap- 
proximately nfa;(l')10~^'^^^'''~^-', while for fields 
smaller than 1 deg^ the sample variance is at least 
nfw(l')10"^-^^(^"^l Table 4 lists uncertainties for 
(p* and Jb derived from subsamples and the angu- 
lar clustering, along with the w(l') and 7 values 
from M. J. I. Brown et al. (in preparation). Sub- 
samples underestimate the uncertainties of ip* and 
js by a factor of ~ 3 at low redshift, so through- 
out the remainder of the paper we use angular 
clustering uncertainties for (f>* and Jb- 

Our 1/Vmax luminosity function values and 
Schechter function fit parameters are provided in 



Tables 5 and 6 respectively. For comparison to the 
prior literature, fits with a ~ —0.5 are also pro- 
vided in Table 6 though for the remainder of the 
paper we use fits where a is a free parameter. As 
shown in Figure 11, Schechter functions provide 
a good approximation of the observed luminosity 
function. 

It is evident from the fits shown in Figure 11 
that the red galaxy luminosity function evolves 
with redshift. In particular, the bright end of the 
luminosity function (i.e., at L > L*) steadily fades 
by ~ 0.9 magnitudes per unit redshift from z — 0.9 
to z = 0.3. While a may decrease slightly with in- 
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Fig. 11. — Luminosity functions for our red galaxy sample. The solid lines are maximum likelihood Schechter 
function fits to the data while the symbols are 1/Vmax estimates of the luminosity function. We use the same 
symbols for our four redshift bins throughout the remainder of the paper. For clarity, we do not plot the la 
uncertainties here but list them in Table 5. For most of our redshift and luminosity bins, these uncertainties 
are on the order of 10%. Within our sample, red galaxies brighter than L* are clearly evolving, while a and 
(p* do not show a strong trend with increasing redshift. 



creasing redshift, the fits computed for each red- 
shift are constrained by different ranges of galaxy 
luminosity and the correlation of a with redshift 
disappears when we only fit the luminosity func- 
tion of Mb — 5\ogh < —19 — 0.92; red galax- 
ies. From Figure 11 and Table 5 it is unclear 
if the luminosity function is or is not undergo- 
ing density evolution. While an apparent decline 
in the space density is measured (see Table 5), 
the 1/Vmax luminosity function bins are correlated 
with each other and Figure 11 does not include 
the uncertainties from cosmic variance, which are 
~ 0.1(/)(M). The normalization of the luminosity 
function, ip* , also does not show a consistent trend 
with redshift. We note, however, that cf)* is very 
sensitive to the measured and assumed values of 
both M* and a. 

In §5.1 and §5.2 we measure the assembly and 
evolution of red galaxies using the luminosity den- 
sity and the luminosity evolution of galaxies at a 
fixed space density threshold. These parameteri- 
zations of the evolving luminosity function are not 
highly correlated with Schechter function param- 



eters such as a. Also, these parameters can be 
easily compared with simple models of red galaxy 
assembly and stellar population evolution, thus 
simplifying the interpretation of the evolving lu- 
minosity function of red galaxies. 

5.1. THE LUMINOSITY DENSITY AND 
THE EVOLVING STELLAR MASS 
WITHIN RED GALAXIES 

We use the luminosity density to measure the 
rest-frame _B-band light emitted by ensembles of 
red galaxies. The luminosity density is strongly 
correlated with the total stellar mass contained 
within all red galaxies, though it is also sensitive 
to the evolving stellar populations of these galax- 
ies. To evaluate the luminosity density for all lu- 
minosities, we integrate over the best-fit Schechter 
functions, so 

Jb - rBL*BT{a + 2). (6) 

We present the luminosity density of all red galax- 
ies as function of redshift in Figure 12. Lumi- 
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Fig. 12. — The B-band luminosity density of red galaxies with redshift. While a Bruzual & Chariot (2003) 
r = 0.6 Gyr model can approximate the colors of red galaxies, it overestimates the evolution of the luminosity 
density. A fit to our data and the 2dFGRS is shown with the dashed blue line. This fit indicates that the 
luminosity density has increased by 36 ± 13% from 2: = to z = 1, rather than the 213% predicted by the 
r — 0.6 Gyr model. If red galaxy stellar populations have faded by ~ 1.24 B-band magnitudes since z — 1, 
the stellar mass contained within the ensemble of red galaxies has approximately doubled over the same 
period. 



nosity density values are listed in Table 7, includ- 
ing values determined with galaxies brighter than 
evolving absolute magnitude thresholds which are 
a function of A/*. For our range of a values, less 
than 25% of the luminosity density is contributed 
by galaxies fainter than our magnitude limits and 
approximately 70% of the luminosity density is 
contributed by galaxies within a magnitude of M* . 
Our measurements are broadly consistent with the 
prior literature, except at z — 0.9 where we find 
a higher luminosity density than COMBO-17 and 
DEEP2. 

To model the evolution of the luminosity den- 
sity, we assume that the stellar populations of 
red galaxies can be approximated by a Bruzual & 
Chariot (2003) r = 0.6 Gyr model with a forma- 
tion redshift of z = 4. This model approximates 
the evolving rest-frame colors of My — Slogft. ~ 
—21 red galaxies, which we plot in Figure 4. The 
T = 0.6 Gyr model fades by 1.24 i?-band mag- 
nitudes between z = I and z — 0, which is con- 
sistent with the observed luminosity evolution of 



both the size-luminosity (Figure 5) and fundamen- 
tal plane relations (e.g., van Dokkum & Stanford 

2003; Treu et al. 2005) of My - 5logh 21 red 

galaxies. We caution that less luminous red galax- 
ies are bluer than the r = 0.6 Gyr model and their 
stellar populations may exhibit more rapid lumi- 
nosity evolution at z < 1 (e.g., Treu et al. 2005; 
van der Wei et al. 2005). The r = 0.6 Gyr model, 
normalized to the 2dFGRS, is plotted in Figure 12 
and grossly overestimates the evolution of the lu- 
minosity density of red galaxies. 

The relationship between log(js) and redshift 
in Figure 12 can be approximated by a straight 
line. If we fit to the Bootes data alone, we find 

log[jB(10^/i i© Mpc"^)] = 7.82(±0.06) 

-K0.22(±0.09) X z, 
(7) 

with a reduced of only 0.2. While the t = 
0.6 Gyr model predicts a 213% increase in the B- 
band luminosity density of red galaxies between 
z = and z = 1, using the Bootes data alone we 
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observe an increase of only 65 ± 34%. 

Although we have the largest 0.2 < z < 1.0 
red galaxy sample currently available, our mea- 
surement of Jb evolution is not sufficiently precise 
to clearly distinguish between significantly difi^er- 
ent descriptions of luminosity function evolution: 
pure luminosity evolution (where Jb evolves in the 
same manner as L*) or luminosity evolution with 
up to 50% density evolution between z = and 
z = 1. While combining data sets can introduce 
systematic errors, the combined data set allows 
the investigation of Jb over a longer redshift base- 
line and reduces the random uncertainties at low 
redshift. The 2dFGRS red galaxy sample is se- 
lected using a criterion based upon principal com- 
ponent analysis of 2dF spectra (Madgwick et al. 
2002), which is strongly correlated with galaxy 
rest-frame color. By fitting Bruzual & Chariot 
(2003) T models to galaxies with SDSS photom- 
etry and 2dFGRS spectroscopy, we find the 2dF 
red galaxy selection criterion corresponds to rest- 
frame U — V > 1.1, which is very similar to our 
criterion for L* red galaxies at low redshift. Using 
the combined data set, we find 

log[iB(10^/i Lq Mpc-3)] = 7.88(±0.02) 

+0.13(±0.04) X z, 
(8) 

which has a reduced of 0.5. Our best fit to the 
combined data sot. shown in Figure 12, yields a 
36 ± 13% increase in the luminosity density of red 
galaxies between z = and z = 1. 

We can infer the rate of stellar mass growth 
within the red galaxy population by comparing the 
evolution of Jb with the luminosity evolution pre- 
dicted by a stellar population synthesis model. We 
assume red galaxy stellar populations have faded 
by 1.24 i?-band magnitudes since 2 = 1, as does 
the T = 0.6 Gyr model. We caution that if we un- 
derestimate the himinosity evolution of the stel- 
lar populations, we underestimate the growth of 
stellar mass within red galaxies, and vice versa. 
If red galaxy stellar populations fade by 1.24 B- 
band magnitudes per unit redshift, then using the 
Bootes data alone we find the stellar mass con- 
tained within the red galaxy population has in- 
creased by 91 ± 39% since z = 1. Using both 
Bootes and the 2dFGRS, we find the stellar mass 
contained within the red galaxy population has 



increased by 131 ± 22% since z = 1. Unless wc 
have grossly overestimated the luminosity evolu- 
tion of red galaxy stellar populations, the stellar 
mass contained within the red galaxy population 
has increased by order unity since z = 1. 

The mild evolution of red galaxy luminosity 
density rules out several models of red galaxy 
evolution. Stellar population synthesis models 
with roughly constant star formation histories ex- 
hibit modest evolution of their i3-band luminosi- 
ties at 2; < 1. However, such star formation his- 
tories result in galaxies with rest-frame colors of 
U — V < 0.5, which is bluer than any of our red 
galaxies. While red galaxy mergers occur at ^ < 1 
(e.g., Lauer 1988; van Dokkum 2005), unless they 
are accompanied by star formation they will only 
redistribute stellar mass already within the red 
galaxy population. Stellar mass must be added to 
the red galaxy population from an outside source, 
namely blue galaxies. 

Blue galaxies must be adding mass to the red 
population due to a decline in their star formation 
rates. Stellar population synthesis models with 
1.5 Gyr < r < 2.5 Gyr move across our rest-frame 
U — V selection criterion for L* galaxies between 
z = 1 and z = 0. Stellar population synthesis 
models can also move across our selection criterion 
if star formation is rapidly truncated, possibly af- 
ter AGN feedback or a merger triggered starburst. 
While a broad range of star formation histories can 
transform blue galaxies into red galaxies, the tight 
color-magnitude relation of red galaxies indicates 
that red galaxies of a given mass and metallicity 
have similar low rates of star formation. These 
galaxies have similar star formation histories over 
many Gyr, or they converge towards similar yet 
low star formation rates upon entering the red 
galaxy population. We will explore some of these 
possibilities in detail in future papers. 

5.2. THE EVOLUTION OF VERY LU- 
MINOUS RED GALAXIES 

Fading star-forming galaxies cannot explain the 
evolution of 4L* red galaxies at z < 1, as star- 
forming galaxies with comparable masses are ex- 
ceptionally rare at these redshifts (e.g.. Bell et al. 
2004). If the stellar mass contained within these 
very luminous red galaxies is evolving at 2 < 1, 
this is presumably due to galaxy mergers. 
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The luminosity density is a poor measure of the 
evolution of 4L* red galaxies. As these galaxies are 
on the exponential part of the luminosity function, 
the measured luminosity density and space den- 
sity of galaxies brighter than an absolute magni- 
tude threshold will have an extremely strong de- 
pendence upon that threshold. For example, in 
Table 7 we hst values of Jb for Mb < M* - 1.5 
galaxies, but M* is correlated with a so the Jb val- 
ues determined with floating and fixed values of a 
differ by up to 35%. Similarly, if the measured lu- 
minosities or the assumed luminosity evolution of 
4L* red galaxies are in error 0.1 magnitudes, the 
measured space density evolution of red galaxies 
can be in error by 30%. If the stellar populations of 
today's L* and AL* galaxies did not evolve in the 
same manner, luminosity density values derived 
with Mb < M* — 1.5 red galaxies can provide a 
misleading picture of 4L* red galaxy evolution. 

We avoid these issues by measuring the evolv- 
ing luminosity function at a fixed space den- 
sity threshold. We parameterize the luminosity 
evolution with M(10~^-^), the absolute magni- 
tude where the space density of red galaxies is 
10~^-^ Mpc^'' mag^^. The most luminous 
red galaxies largely lie along the color magnitude 
relation, so M(10~^'^) is insensitive to details of 
red galaxy selection criteria and it is thus easy to 
compare the results of various surveys. Schechter 
functions provide a good fit to the red galaxy lu- 
minosity fmiction when the space density is above 
10-'*-^ h-^ Mpc"^ mag-^ (e.g., Madgwick et al. 
2002), so we use best-fit Schechter functions to 
derive M(10~^'^) from our data and the literature 
(when available). Our estimates of M(10~^'^) 
are provided in Table 6, and plotted in Figure 13 
along with values derived from the literature. We 
find Mb (10"^ ''') steadily increases with increasing 
redshift. A straight line fit to the Bootes data 
alone has a reduced of 1.1 and A'/b(10^"^'^) 
brightens by 0.87 ± 0.06 i?-band magnitudes from 
2; = to 2; = 1. While there are some discrep- 
ancies between the surveys at z ~ 0.9. whic:li we 
discuss in §6, there is broad agreement between 
our results and the literature at 2; < 0.8. If we 
fit to both the Bootes and 2dFGRS data, we 
find M(10-3-5) brightens by 0.87 ± 0.05 S-band 
magnitudes between 2: = and 2; = 1. 

We use two simple models to characterize the 
evolution of Mb(10~'^-^). Both models are in- 



tended to be illustrative rather than precise de- 
scriptions of galaxy evolution. Our simplest model 
has no galaxy mergers and assumes that the star 
formation history of all red galaxy progenitors is 
described by a Bruzual & Chariot (2003) t = 
0.6 Gyr model. This r model has negligible star 
formation at z < 1 and approximates the evolu- 
tion of the color-magnitude and size-luminosity re- 
lations of ~ AL* red galaxies (Figures 4 and 5). 
As mergers of red galaxies do occur at z < 1 (e.g., 
Lauer 1988; van Dokkum 2005), this model pro- 
vides an upper bound for the fading of Mb(10~^'^) 
from 2 = 1 to 2; = 0. 

Our second model of Mb(10~^'^) assumes 
galaxy mergers can be described by the Sheth 
& Tormen (1999) halo mass function, which is 
similar to the formalism of Press & Schechter 
(1974). We assume one galaxy per dark matter 
halo, a constant ratio of baryonic matter to dark 
matter in each halo, and a r = 0.6 Gyr model 
star formation history. From the halo mass func- 
tion we can determine the space density of halos 
more massive than nimin, and we know the space 
density of galaxies brighter than Mb {10^^-^) is 
~ 10~^/i^Mpc~^. We therefore assume that the 
growth of baryons within Mb(10^'^'^) galaxies is 
proportional to the evolving mmin value corre- 
sponding to a halo space density of lO-^/i^Mpc"^. 
We find rrimin > IO^^Mq at z < 1, which is com- 
parable to the masses of galaxy groups. This sim- 
ple model is clearly flawed, as group and cluster 
halos contain multiple galaxies, and the infall of 
galaxies into these halos may not result in galaxy 
mergers. As this model overestimates the rate of 
galaxy mergers, it should provide a lower bound 
for the evolution of Mb(10-^-^). 

Our MBilO^"^'^) models are plotted in Fig- 
ure 13, along with Mb{10^'^'^) values derived from 
this work and the literature. The model with- 
out mergers clearly provides the better approx- 
imation to the data. This model is offset from 
the data by only 0.21 magnitudes a,t z — 0.7 and 
0.31 magnitudes at z = 0.9. If the stellar pop- 
ulations of red galaxies fade in the same manner 
as the T — 0.6 Gyr model, ~ 75% of the stellar 
mass contained within today's AL* galaxies was 
already in place within these galaxies by 2; = 0.9. 
If red galaxy stellar populations fade by less than 
1.24 B-band magnitudes from 2: = 1 to 2 = 0, an 
even higher percentage of the stellar mass of AL* 
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Fig. 13. — The luminosity evolution of ~ 4L* red galaxies. We have parameterized the evolution with 
MsilO^"^'^), the absolute magnitude where the space density of red galaxies is 10~'^'^ Mpc"'' mag""'^. 
Random uncertainties for the Bootes data are smaller than the size of the data points. Simple models with 
and without stellar mass growth via galaxy mergers are plotted, and these models are described in detail in 
§5.2. Both models are normalized to the 2dFGRS and assume the stellar populations of red galaxies fade 
by 1.24 _B-band magnitudes per unit redshift. Our simple merging model, based upon the Sheth & Tormen 
(1999) mass function, underestimates the space density of 4L* red galaxies at z > 0.5. While our model 
without mergers does not fit the data, it is only offset from our data by 0.21 magnitudes at z = 0.7, and we 
therefore conclude that ~ 80% of the stellar mass contained within today's ~ 4L* red galaxies was already 
in place at z = 0.7. 



red galaxies is already in place prior to z = 0.9, 
and vice versa. While we measure a higher space 
density of very luminous red galaxies at z = 0.9 
than the literature, there is broad agreement at 
lower redshifts. Even if the Bootes measurements 
at z = 0.9 are in error, the space density of 4L* 
red galaxies from the literature is approximated 
by our model at z < 0.7. Roughly 80% of the 
stellar mass within today's AL* red galaxies must 
already be in place by z = 0.7. 

We can see in Figure 13 that the simple merg- 
ing model, based upon the Sheth & Tormen (1999) 
halo mass function for a ACDM cosmology, dif- 
fers significantly from the observed evolution of 
M(10~'^'^). The simplifying assumptions used in 
this model may explain the offset from the data. 
We assumed one galaxy per halo, while the most 
luminous red galaxies often reside within galaxy 
groups and clusters. An accurate model of these 



galaxies requires a description of how these galax- 
ies occupy subhalos within groups and clusters, 
and how these subhalos merge over cosmic time. 
Such a model must also have negligible star for- 
mation at z < 1, and clearly there is a physical 
process that is extremely effective at truncating 
star formation in 4L* red galaxies. Issues such 
as the halo occupation distribution function (e.g., 
Berlind & Weinberg 2002), gas heating (e.g., Naab 
et al. 2005), and AGN feedback (e.g.. Bower et al. 
2006; Croton et al. 2006; Hopkins et al. 2006) are 
clearly important for understanding galaxy evolu- 
tion and need to be studied in more detail. 

6. COMPARISON TO OTHER RED 
GALAXY LUMINOSITY FUNCTIONS 

There have been a variety of conclusions about 
red galaxy evolution over the past decade (see 
Faber et al. 2005), and some of this can be ex- 
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Fig. 14. — A comparison of our luminosity functions with those from the 2dFGRS and SDSS (Madgwick at al. 
2002; Blanton 2005). To plot the 2dFGRS and SDSS luminosity functions, we have adopted Bj - B = 0.15 
and g^B — ^Voga- There is a steady brightening of > L* red galaxies with increasing redshift. 
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Fig. 15. — A comparison of our luminosity functions with those of COMBO-17 (Faber et al. 2005), DEEP2 
(Willmer et al. 2005; Faber et al. 2005), and the reanalysis of DEEP2 by Blanton (2005). For clarity we 
don't plot the VVDS luminosity function of Zucca et al. (2005), as their non-evolving red galaxy selection 
criterion differs greatly from the others displayed in this plot. The luminosity functions broadly agree at 
z — 0.5, but at z = 0.9 there are significant disagreements between each of the surveys, particularly for 
J > 21 red galaxies. Potential causes for these discrepancies are discussed in §6. 
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plained by the different behavior of L* and 4L* 
red galaxies (§5.1 and §5.2 respectively). Bundy 
et al. (2005) also find that the evolution of red 
galaxies is a function of stellar mass, and this has 
been confirmed by recent analyzes of COMBO-17 
stellar mass and luminosity functions (Borch et al. 
2006; Cimatti, Daddi, & Renzini 2006). The lumi- 
nosity density of red galaxies is heavily weighted 
towards L* galaxies and evolves slowly at z < 1. 
After accounting for fading stellar populations, we 
and much of the recent literature (e.g.. Bell et al. 
2004; Faber et al. 2005) conclude that the stellar 
mass contained within the ensemble of red galaxies 
has steadily increased since z = 1. In contrast, the 
evolution of AL* red galaxies differs only slightly 
from pure luminosity evolution. However, unlike 
some recent studies (e.g., Bundy et al. 2005; Borch 
et al. 2006; Cimatti et al. 2006), we do see evi- 
dence for the ongoing assembly of the most mas- 
sive red galaxies, albeit at a rate that produces 
little growth of AL* red galaxy stellar masses at 
z < 1. As red galaxy evolution is a function of lu- 
minosity, it is not surprising that various authors 
using different techniques have reached a variety of 
conclusions, even when using the same galaxy sur- 
veys (e.g.. Bell et al. 2004; Willmer et al. 2005; 
Faber et al. 2005; Blanton 2005; Bundy et al. 
2005; Borch et al. 2006; Cimatti et al. 2006). 

Figures 14 and 15 show our evolving luminosity 
functions are broadly consistent with the 2dFGRS, 
SDSS, COMBO-17 and DEEP2 at z < 0.8. There 
is a steady increase of the bright end of the lumi- 
nosity function with redshift. If we combine our 
results with those of the 2dFGRS and SDSS, there 
is a gradual decline of (/)* with increasing redshift. 
However, we show in Figure 15 that we measure a 
higher space density of z = 0.9 red galaxies than 
DEEP2 or COMBO-17. Since these studies have 
had a significant impact upon the field, it is impor- 
tant to understand why these differences occur. 

An obvious difference between our work and 
some of the recent literature is our use of photo- 
metric redshifts. Photometric redshift errors can 
scatter numerous low luminosity objects into high 
luminosity bins, and thus alter the shape of the lu- 
minosity function (e.g., Brown, Webster, & Boyle 
2001; Chen et al. 2003). To account for this we 
convolved our evolving luminosity function by the 
expected luminosity and volume errors resulting 
from our photometric redshifts. We use the pho- 
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Fig. 16. — The effect of photometric redshift errors 
on the measured z = 0.9 luminosity function. We 
use the photometric redshift uncertainties listed in 
Table 1, except when there are ten or less spectro- 
scopic redshifts, in which case we extrapolate from 
brighter magnitudes by assuming the uncertainties 
increase by 50% per unit magnitude. The effect of 
photometric redshift errors upon the luminosity 
function is only significant for M^ — 51ogft- < —22.5 
galaxies, where the space density is overestimated. 
In Figure 11 the excess of very luminous galaxies 
compared to a Schechter function may be caused 
by photometric redshift errors. Even if our photo- 
metric redshift uncertainties are 50% larger those 
listed in Table 1, the impact on our results and 
conclusions is small as the measured luminosity 
function of -21.8 < Mb - 5logh < -19.5 red 
galaxies remains largely unchanged. 

tometric redshift errors listed in Table 1, except 
at faint magnitudes with ten or less spectroscopic 
redshifts, where we extrapolate from brighter mag- 
nitudes by assuming the errors increase by 50% per 
unit magnitude. In Figure 16 the largest error is 
a ~ 0.2 magnitude offset for rare lOL* galaxies at 
z = 0.9. In Figure 11, there are slightly more lOL* 
galaxies than expected from our best-fit z = 0.9 
Schechter function, but these are a tiny fraction 
of our sample and do not affect our results, which 
rely upon less luminous galaxies. 

Compared to the other luminosity functions in 
Table 6, the VIRMOS-VLT Deep Survey (VVDS; 
Zucca et al. 2005) find a lower space density 
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of high rcdshift and low himinosity red galax- 
ies. As discussed by Zucca et al. (2005), this 
is not unexpected as the WDS selection criteria 
do not model the evolution and tilt of the color- 
magnitude relation. Wolf et al. (2003) also ob- 
served a rapid decline of (f)* with redshift for a 
sample of red galaxies selected from COMBO-17 
with a non-evolving criterion. Similar selection 
effects may occur in other surveys by accident at 
z > 1, where the evolving color- magnitude relation 
is relatively difficult to measure and red galaxy se- 
lection criteria could unintentionally intercept the 
color-magnitude relation. 

Both the WDS and COMBO-17 use SExtrac- 
tor MAG_AUTO to determine total magnitudes, 
and accuracy of the SDSS DRl and 2dFGRS pho- 
tometry has been verified with deep CCD pho- 
tometry making use of MAG_AUTO (Cross et al. 
2004). MAG_AUTO is known to underestimate 
the luminosities of galaxies with de Vaucouleurs 
(1948) profiles by 0.1 magnitudes when the half- 
light radius is several times larger than the stellar 
full width at half maximum (Cross et al. 2004) . It 
is therefore plausible that the brightest red galax- 
ies in the 2dFGRS have had their luminosities un- 
derestimated by ~ 10%. As illustrated in Figure 1, 
MAG_AUTO can also have systematics of several 
tenths of a magnitude at / > 21, even when the 
catalogs are 85% complete. We therefore specu- 
late that the WDS and COMBO-17 may have 
systematic errors on the order of 0.2 magnitudes 
at the faint end of their luminosity functions. 

While our z ~ 0.9 luminosity functions and 
those derived from DEEP2 (Willmer et al. 2005; 
Fabcr et al. 2005: Blanton 2005) arc broadly sim- 
ilar in shape and normalization, there is a sys- 
tematic offset of ~ 0.2 magnitudes. This offset 
is present even at high luminosities, which are 
largely insensitive to the selection effects discussed 
elsewhere in this paper. As discussed in §3.4, we 
correct for flux outside of our aperture using a 
model derived from the size-luminosity relation. 
Such corrections are typically not applied to sam- 
ples in the literature. 

We have estimated the fraction of galaxy flux 
that falls beyond the photometric aperture used 
by the DEEP2 survey. DEEP2 photometry uses 
the larger of a 2" diameter aperture or a 6a aper- 
ture, where a is defined by a Gaussian fit to the 
object profile (Coil et al. 2004). DEEP2 pho- 




I— band magnitude (Vega) 

Fig. 17. — The DEEP2 aperture diameter for 
R - I > 1.2, 0.8 < 2: < 1.0 galaxies from the 
first DEEP2 data release. The apparent color cut 
effectively selects galaxies with red rest-frame col- 
ors in this redshift range. On the right we list the 
expected magnitude offsets for galaxies with a 0.5" 
half-light radius observed in 0.8" seeing and mea- 
sured with the DEEP2 aperture. As the DEEP2 
photometry is zero-pointed with stars measured in 
a 2" aperture, larger apertures overestimate the 
luminosities of point sources and very large aper- 
tures can overestimate galaxy luminosities. How- 
ever, wo find that DEEP2 photometry typically 
underestimates the total luminosities of z ~ 0.9 
red galaxies. 



tometry is zero-pointed with SDSS photometry of 
stars, where the SDSS uses PSF magnitudes while 
DEEP2 uses apertures of ~ 2" diameter. Because 
of this, a point source measured with a > 2" diam- 
eter aperture will have its flux systematically over- 
estimated. Similarly, the photometry of galaxies 
will have systematic errors which are a function of 
half-light radius and aperture size. In Figure 17, 
we plot the DEEP2 aperture diameter as a func- 
tion of apparent magnitude for 0.8 < z < 1.0 red 
galaxies selected from the first DEEP2 data re- 
lease. At these redshifts, an / ^ 21.5 galaxy has 
a half-light radius of ~ 0.5" (~ 3/i~^kpc), and we 
list the systematic error for such a galaxy mea- 
sured with the DEEP2 aperture on the right of 
Figure 17. Red galaxies at 2; ~ 0.9 may have 
their luminosities systematically underestimated 
by ~ 0.15 magnitudes by DEEP2, and we there- 
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Fig. 18. — The difference between absolute magnitudes and rest-frame colors determined with Bruzual & 
Chariot (2003) models and Kinney et al. (1996) templates. For this plot we have only included Ms — Slog/i < 
— 19.5 red galaxies. To produce a continuous model of the galaxy locus, we interpolated between the Kinney 
et al. (1996) SEDs of elliptical, SO, Sa, Sb, and Sc galaxies. The red solid hne shows the mean offset as a 
function of redshift. While the i3-band absolute magnitudes are relatively insensitive to the choice of SED 
models or templates, the rest-frame U — V colors show systematics on the order of 0.2 magnitudes. If we 
have overestimated the evolution of U — V with redshift, our principal conclusions remain unchanged. 



fore conclude that different measures of total mag- 
nitudes contribute to the difference between our 
work and the prior literature. 

Different surveys use different SED models and 
templates to derive rest-frame galaxy properties 
from observed photometry. Figure 18 compares 
rest-frame magnitudes and colors derived with the 
Bruzual & Chariot (2003) r models and Kinney 
et al. (1996) templates for our red galaxy sam- 
ple. The absolute magnitudes exhibit small dif- 
ferences, as the observed /-band is very close to 
the rest-frame i3-band at z = 0.9. In contrast, 
the rest-frame U — V colors show systematics as 
large as 0.2 magnitudes. This will have minimal 
impact on the bright end of the luminosity func- 
tion, as the most luminous red galaxies lie along 
the color-magnitude relation, which can be iden- 
tified empirically. The greatest impact from SED 
errors will be at lower luminosities, as SED errors 
will shift red galaxy selection criterion relative to 
the the color- magnitude relation. A 0.1 magnitude 
shift in our U ~V selection criterion results in our 
(j)* and Jb values changing by ^ 25%. 

If the Kinney et al. (1996) templates are bet- 



ter models of z ^ 0.9 red galaxy SEDs than the 
Bruzual & Chariot (2003) t models, our principal 
conclusions remain unchanged. The U ^ V color- 
magnitude relation still evolves, but the mean 
color and luminosity of red galaxy stellar popu- 
lations evolves less rapidly than predicted by the 
T = 0.6 Gyr model. Even if our z = 0.9 js mea- 
surement has a 25% error, the stellar mass con- 
tained within the red galaxy population must in- 
crease with decreasing redshift. If we reduce the 
luminosity evolution of red galaxy stellar popula- 
tions, a revised Mb{W~^'^) model without merg- 
ers will provide a better fit to the observations in 
Figure 13 than our current model. 

7. SUMMARY 

We have measured the evolution and assembly 
of red galaxies over the last 8 Gyr, using the evolv- 
ing i3-band luminosity function of 0.2 < z < 1.0 
red galaxies. Our sample of 39599 0.2 < z < 1.0 
red galaxies, selected from 6.96 deg^ of imaging, 
is an order of magnitude larger in size and vol- 
ume than comparable z ~ 0.9 red galaxy sam- 
ples. Accurate photometric redshifts were deter- 
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mined using the empirical ANNz photometric red- 
shift code and photometry from both the NOAO 
Deep Wide-Field and Spitzer IRAC shallow sur- 
veys. The accuracy of these redshifts has been 
verified with spectroscopy, and is comparable to 
the best broad-band photometric redshifts in the 
literature. 

We find the color-magnitude and size- luminosity 
relations of red galaxies evolve with redshift, and 
this evolution can be approximated by a Bruzual 
& Chariot (2003) r = 0.6 Gyr stellar population 
synthesis model. The size-luminosity relation pre- 
dicts a significant fraction of the galaxy flux is at 
radii of several arcseconds, even at ^; ~ 0.9. Lu- 
minosity functions which do not account for this 
underestimate the space density of red galaxies at 
high redshift, and overestimate the assembly rate 
of red galaxies at ^; < 1. 

We find the luminosity density, js, of red galax- 
ies increases by 36 ± 13% from z = to 2 = 1. In 
contrast, the i?-band luminosity of a r = 0.6 Gyr 
stellar population model increases by 213% over 
the same redshift range. If red galaxy stellar pop- 
ulations fade by 1.24 magnitudes from 2; = 1 to 
z = 0, the stellar mass contained within red galax- 
ies has approximately doubled over the past 8 Gyr. 
Blue galaxies are being transformed into < L* red 
galaxies at 2 < 1, after a steady decline or rapid 
truncation of their star formation. 

The evolution of the most luminous red galax- 
ies at z < 1 cannot be dominated by fading 
blue galaxies, as blue galaxies with large stel- 
lar masses are exceptionally rare at 2; < 1. 
We have measured the evolution of AL* red 
galaxies using Mb(10~^'^), the absolute magni- 
tude corresponding to a fixed space density of 
10~^'^ Mpc~^ mag~^. Our measurements of 
AfB(10"^-''), along with those we derived from 
the literature, show that the stellar masses of the 
most luminous red galaxies evolve slowly at 2; < 1. 
A Bruzual & Chariot (2003) t = 0.6 Gyr model, 
with little stellar mass evolution and normalized 
to the 2dFGRS, only overestimates Mb{W-^-^) 
by ~ 0.21 magnitudes at z = 0.7. We therefore 
conclude that ~ 80% of the stellar mass contained 
within today's Mg — 51og/i = — 21 red galaxies was 
already in place within these galaxies at z = 0.7. 
While red galaxy mergers have been reported in 
the prior literature, such mergers do not rapidly 
increase the stellar masses of 4Z/* red galaxies 



between z = 1 and the present day. 
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Table 1 

Measured Photometric Redsiiift I.Tncertaixties For Red Galaxies 





0.2 < z < 0.4 


0.4 < 2 < 0.6 


0.6 < z < O.i 


i 


O.i 


i < 2 < 1.0 


4" Apparent 


Mb - 


51og/» < 


-17.5 


Mb - 


51og/i < 


-18.0 


Mb - 


5logh < - 


-19.0 


Mb - 


51og/i < - 


19.5 


Magnitude 


68.7% 


90% 


N 


68.7% 


90% 


N 


68.7% 


90% 


N 


68.7% 


90% 


N 


17.0 < / < 18.0 


0.033 


0.060 


(170) 




















18.0 < / < 19.0 


0.032 


0.063 


(1110) 


0.026 


0.051 


(132) 


0.025 


0.025 


(1) 








19.0 < / < 20.0 


0.039 


0.074 


(1007) 


0.032 


0.058 


(949) 


0.029 


0.053 


(192) 


0.027 


0.179 


(10) 


20.0 < / < 21.0 


0.049 


0.088 


(166) 


0.036 


0.066 


(340) 


0.040 


0.072 


(240) 


0.058 


0.096 


(50) 


21.0 < I< 22.0 


0.042 


0.363 


(4) 


0.229 


3.066 


(5) 


0.103 


0.131 


(14) 


0.045 


0.090 


(28) 


22.0 < / < 23.5 














0.162 


0.271 


(5) 


0.127 


0.254 


(7) 



Table 2 

Color cuts used to exclude stars, quasars and z > 1 galaxies from the sample. 



Magnitude and/or color range 


Cut 




Bw < 26.5 


Bw - R< 1.00 




Biv < 26.5 


Bw - R> 4.20 




Bw > 26.5 


26.5 - R> 4.20 




Bw < 26.5 


Bw - R> 2.55 + 2.75 X {R - I - 0.55) 




Bw > 26.5 


26.5 - R> 2.55 + 2.75 x (R - I - 0.55) 




/ < 23.5 


R- I < 0.55 




/ < 23.5 


R- I > 1.70 




Bw < 26.5 


Bw - R < 2.0 - {R - I) 




Bw < 26.5 


Bw - R < 5.0 X (R - I - 1.2) 




[3.6] < 19.5 


/ - [3.6] < 2.2 




[3.6] > 19.5 


/ - 19.5 < 2.2 




[3.6] < 19.5 and 0.75 < R - I < 1.05 


R- I > 0.75 + 0.375 x (/ - [3.6] - 2.2) 




[3.6] > 19.5 and 0.75 < ii - / < 1.05 


R- I > 0.75 + 0.375 x (/ - 19.5 - 2.2) 




[3.6] < 19.5 and 1.05 < R - / < 1.70 


R- I > 1.05 + 0.60 X (/ - [3.6] - 3.0) 




[3.6] > 19.5 and 1.05 < B. - I < 1.70 


R- I > 1.05 + 0.60 X (/ - 19.5 - 3.0) 




[3.6] < 19.5 and [;!.(>] < 18.5 


[3.6] - [4.5] > 0.6 




[3.6] < 19.5 and [3.6] < 18.5 


[3.6] - [4.5] < -0.6 






Table 3 






Red Galaxy Number Counts 




4" Apparent 0.2 < z < 0.4 


0.4 < z < 0.6 0.6 < 2 < 0.8 


0.8 < z < 1.0 


Magnitude Mb - 51og/t < -17.5 Mb - 5logh < -18.0 Mb - 51og/t < -19.0 


Mb - 5loeh < -19.5 


17.0 < / < 18.0 189 






18.0 < / < 19.0 1272 


141 1 




19.0 < / < 20.0 2365 


2131 400 


14 


20.0 < / < 21.0 1719 


4591 3320 


977 


21.0 < I < 22.0 438 


3512 5738 


5466 


22.0 < I < 23.5 


608 1359 


5358 



Table 4 

Comparison of subsample and clustering uncertainties for 6* and jp, 



z range 


0,(1') 


7 


(/>* X 

Best fit 
Value 


10^(/i-'Mpc"-^ 
Subsample 
Icr 


mag ^) 
Clustering 

1(7 


Best fit 
Value 


(lO^/i Lq Mpc 
Subsample 
lu 


Clustering 
Iff 


0.20 < z < 0.40 


0.74 ± 0.05 


1.94 ± 0.05 


8.45 


±0.38 


±1.07 


8.1 


±0.3 


±1.0 


0.40 < z < 0.60 


0.51 ± 0.03 


2.02 ± 0.07 


7.61 


±0.34 


±0.68 


8.5 


±0.4 


±0.8 


0.60 < z < 0.80 


0.28 ± 0.02 


2.06 ± 0.06 


5.71 


±0.33 


±0.35 


9.2 


±0.4 


±0.6 


0.80 < z < 1.00 


0.24 ±0.02 


1.93 ±0.05 


6.35 


±0.21 


±0.47 


10.7 


±0.3 


±0.8 



27 



Table 5 

Red Galaxy 1/Vmax Luminosity Function with subsample uncertainties 



Absolute 




Luminosity Function {h' 


' Mpc~^ mag~^) 




Magnitude 


0.2 < z < 0.4 


OA < z < 0.6 


0.6 < z < 0.8 


0.8 < z < 1.0 



-17.75 


< 


Mb 


— 5Iog/i 


< 


-17 


.50 


1 


.99 ± 0. 


10 


X 


10" 


-3 






























-18.00 


< 


Mb 


— 5Iog/i 


< 


-17 


.75 


1 


.99 ± 0. 


13 


X 


10" 


3 






























-18.25 


< 


Mb 


— 51og/i 


< 


-18 


.00 


2 


.07 ± 0. 


10 


X 


10" 


3 


1 


.94 ± 0.12 


X 10" 


-3 






















-18.50 


< 


Mb 


— 51og/i 


< 


-18 


.25 


2 


.51 ± 0. 


18 


X 


10" 


-3 


2 


.03 ± 0.13 


X 10" 


-3 






















-18.75 


< 


Mb 


— 51og/i 


< 


-18 


.50 


2 


.51 ± 0. 


13 


X 


10" 


-3 


2 


.21 ± 0.11 


X 10" 


-3 






















-19.00 


< 


Mb 


— 51og/i 


< 


-18 


.75 


2 


.98 ± 0. 


17 


X 


10" 


-3 


2 


.71 ± 0.16 


X 10" 


-3 






















-19.25 


< 


Mb 


— 5Iog/i 


< 


-19, 


.00 


3 


.00 ± 0. 


17 


X 


10" 


-3 


2 


.68 ± 0.13 


X 10" 


-3 


2.22 ± 0, 


,11 X 


10" 


-3 














-19.50 


< 


Mb 


— 5Iog/i 


< 


-19, 


.25 


3 


.25 ± 0. 


18 


X 


10" 


-3 


2 


.78 ± 0.11 


X 10" 


3 


2.36 ± 0, 


,11 X 


10" 


-3 














-19.75 


< 


Mb 


— 5Iog/i 


< 


-19, 


.50 


2 


.78 ± 0. 


14 


X 


10" 


-3 


2 


.69 ± 0.10 


X 10" 


3 


2.34 ± 0, 


,11 X 


10" 


-3 


2 


.27 ± 0, 


,08 


X 


10" 


-3 


-20.00 


< 


Mb 


— blogh 


< 


-19 


.75 


2 


.55 ± 0. 


14 


X 


10" 


3 


2 


.56 ± 0.18 


X 10" 


-3 


2.29 ± 0, 


,10 X 


10" 


-3 


2 


.46 ± 0, 


,10 


X 


10" 


-3 


-20.25 


< 


Mb 


— 51og/i 


< 


-20, 


.00 


2 


.17 ± 0. 


14 


X 


10" 


-3 


2 


.03 ± 0.12 


X 10" 


-3 


2.02 ± 0, 


,14 X 


10" 


-3 


2 


.44 ± 0, 


,07 


X 


10" 


-3 


-20.50 


< 


Mb 


— 51og/i 


< 


-20, 


.25 


1 


.64 ± 0, 


,08 


X 


10" 


-3 


1, 


.68 ± 0.07 


X 10" 


-3 


1.80 ± 0, 


,06 X 


10" 


-3 


1, 


.97 ± 0, 


,07 


X 


10" 


-3 


-20.75 


< 


Mb 


— 51og/i 


< 


-20, 


.50 


9, 


.87 ± 0. 


,76 


X 


10" 


-4 


1, 


.24 ± 0.09 


X 10" 


-3 


1.24 ± 0, 


,08 X 


10" 


-3 


1, 


.61 ± 0, 


.07 


X 


10" 


-3 


-21.00 


< 


Mb 


— 51og/i 


< 


-20 


.75 


7, 


.35 ±0. 


96 


X 


10" 


-4 


8, 


.16 ±0.67 


X 10" 


-4 


9.95 ± 0, 


.65 X 


10" 


-4 


1 


.33 ± 0, 


.07 


X 


10" 


-3 


-21.25 


< 


Mb 


— 5\ogh 


< 


-21 


.00 


3 


.37 ± 0. 


,47 


X 


10" 


-4 


5, 


.12 ±0.46 


X 10" 


-4 


6.44 ± 0, 


.34 X 


10" 


-4 


8, 


.97 ±0, 


.46 


X 


10" 


-4 


-21.50 


< 


Mb 


— olog/j. 


< 


-21 


.25 


1 


.96 ± 0. 


22 


X 


10" 


-4 


2 


.62 ±0.37 


X 10" 


-4 


4.66 ±0, 


.41 X 


10" 


-4 


5, 


.94 ± 0, 


.25 


X 


10" 


-4 


-21.75 


< 


Mb 


— 51og/i 


< 


-21 


.50 


3 


.76 ± 2. 


03 


X 


10" 


-5 


1 


.00 ±0.19 


X 10" 


-5 


2.20 ± 0, 


.22 X 


10" 


-4 


3 


.06 ± 0, 


.22 


X 


10" 


-4 


-22.00 


< 


Mb 


— 51og/i 


< 


-21 


.75 


2, 


.38 ± 2. 


,41 


X 


10" 


-5 


4, 


.76 ± 1.09 


X 10" 


-5 


9.47 ± 1, 


.97 X 


10" 


-5 


1 


.34 ± 0, 


.14 


X 


10" 


-4 


-22.25 


< 


Mb 


— 51og/i 


< 


-22 


.00 














1 


.19 ±0.45 


X 10" 


-5 


3.21 ± 0, 


.63 X 


10" 


-5 


5, 


.30 ± 0, 


.74 


X 


10" 


-5 


-22.50 


< 


Mb 


— 51og/i 


< 


-22, 


.25 
















< 1.78 X 


10-^ 




9.16 ± 3, 


.08 X 


10" 


-6 


1 


.96 ± 0, 


.54 


X 


10" 


-5 


-22.75 


< 


Mb 


— 51og/i 


< 


-22, 


.50 
















< 2.20 X 


10-^ 




1.53 ± 1, 


.35 X 


10" 


-6 


4, 


.61 ± 2, 


.42 


X 


10" 


-6 


-23.00 


< 


Mb 


— 51og/i 


< 


-22, 


.75 
















< 3.02 X 


10-^ 




1.53 ± 1, 


.61 X 


10" 


-6 


1 


.15 ± 1, 


.03 


X 


10" 


-6 
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Table 6 

Recent S-band luminosity functions of red galaxies 



Survey* z range Ngaiaxy M* - Slog/t'' Mb (lO"^-^) - 51og/(= (ji*(/(^Mpc-=* mag"^) a'^ 



Bootes 


0.20 < 


z 


< 0, 


.40 


598.3 


-19, 


.54 


± 0. 


,05 


-21.10 ± 0.03 


8.45 ± 1, 


,07 X 10" 


-3 


-0.28 ± 


.04 


Bootes 


0.40 < 


z 


< 0, 


,60 


10983 


-19, 


,72 


± 0, 


,04 


-21.25 ± 0.03 


7.61 ± 0, 


,68 X 10" 


-3 


-0.28 ± 


.05 


Bootes 


0.60 < 


z 


< 0, 


.80 


10817 


-20, 


,16 


± 0, 


,05 


-21.47 ± 0.02 


5.71 ± 0, 


,35 X 10" 


-3 


-0.55 ± 


.06 


Bootes 


0.80 < 


z 


< 1, 


.00 


11816 


-20, 


,21 


± 0, 


,03 


-21.61 ± 0.02 


6.35 ± 0, 


,47 X 10" 


-3 


-0.43 ± 0, 


.05 


Bootes 


0.20 < 


z 


< 0, 


.40 


5983 


-19, 


,78 


±0, 


.02 


-21.19 ±0.03 


7.16 ± 0, 


,90 X 10" 


-3 


-0.50 




Bootes 


0.40 < 


z 


< 0, 


.60 


10983 


-19, 


,92 


±0, 


.02 


-21.30 ±0.03 


6.60 ±0. 


,59 X 10" 


-3 


-0.50 




Bootes 


0.60 < 


z 


< 0, 


.80 


10817 


-20, 


,12 


±0, 


.02 


-21.46 ±0.02 


5.87 ±0. 


,36 X 10" 


-3 


-0.50 




Bootes 


0.80 < 


z 


< 1. 


.00 


11816 


-20, 


,26 


±0, 


.01 


-21.62 ± 0.02 


6.17 ±0. 


,45 X 10" 


-3 


-0.50 





2dFGRS z< 0.15 27540 -19.43 ±0.05 -20.93 9.9 ± 0.5 X 10"^ -0.54 ± 0.02 



COMBO-17 





.2 


< 


z 


< 


0. 


,4 


1096 


-19 


.86 


± 


0. 


,16 


-21, 


,23 


6. 


38 


± 


2. 


,47 


X 


10" 


-3 


-0, 


,50 


COMBO-17 





.4 


< 


z 


< 


0, 


,6 


1179 


-20 


.00 


± 


0. 


,11 


-21, 


,33 


5. 


82 


± 


0, 


,94 


X 


10" 


-3 


-0, 


,50 


COMBO-17 


0, 


.6 


< 


z 


< 


0, 


,8 


1431 


-20, 


.33 


± 


0. 


12 


-21, 


,62 


5. 


17 


± 


0, 


,20 


X 


10" 


-3 


-0, 


,50 


COMBO-17 


0, 


.8 


< 


z 


< 


1, 


,0 


892 


-20 


.41 


± 


0. 


14 


-21, 


,53 


3. 


46 


± 


0, 


,15 


X 


10" 


-3 


-0, 


,50 


COMBO-17 


1 


.0 


< 


z 


< 


1, 


,2 


256 


-20 


.81 


± 


0. 


16 


-21, 


,44 


1. 


55 


± 


0, 


,34 


X 


10" 


-3 


-0, 


,50 


DEEP2 





.2 


< 


z 


< 


0, 


,4 


109 


-20 


.25 


± 


0, 


18 


-21, 


,53 


4. 


97 


± 


0, 


,48 


X 


10" 


-3 


-0, 


,50 


DEEP2 





.4 


< 


z 


< 


0, 


,6 


173 


-20 


.20 


± 


0, 


12 


-21, 


,40 


4. 


13 


± 


0, 


,19 


X 


10" 


-3 


-0, 


,50 


DEEP2 





.6 


< 


z 


< 


0, 


,8 


196 


-20 


.42 


± 


0, 


,06 


-21, 


,60 


3. 


98 


± 


0, 


,31 


X 


10" 


-3 


-0, 


,50 


DEEP2 





.8 


< 


z 


< 


1 


,0 


535 


-20 


.34 


± 


0, 


,05 


-21 


,41 


3. 


13 


± 


0, 


,11 


X 


10" 


-3 


-0, 


,50 


DEEP2 


1 


.0 


< 


z 


< 


1 


,2 


178 


-20 


.67 


± 


0, 


,08 


-21 


,29 


1. 


58 


± 


0, 


,25 


X 


10" 


-3 


-0, 


,50 



VVDS 


0.2 < z < 0.4 


65 


-20.27±«;il 


-21.66 


5.15 ± 0.64 X 10"^ 


-0.29 


VVDS 


0.4 < z < 0.6 


106 


-20.49±°;^ 


-21.67 


3.12 ± 0.30 X 10"^ 


-0.29 


VVDS 


0.6 < z < 0.8 


197 


-20.22±°°» 
-20.73±°ii 


-21.46 


3.53 ± 0.35 X 10-^ 


-0.29 


VVDS 


0.8 < z < 1.0 


164 


-21.76 


2.36 ±0.18 X 10-^ 


-0.29 


VVDS 


1.0 < z < 1.2 


114 


-20.53±l\l 


-21.57 


2.39 ±0.22 X 10-=* 


-0.29 



^Bootes (this work), 2dFGRS (2dF Galaxy Redshift Survey; Madgwick et al. 2002), COMBO-17 (Classifying Objects by Medium- 
Band Observations - a spectrophotometric 17-filter survey; Bell et al. 2004; Faber et al. 2005), DEEP2 (Deep Extragalactic Evolution 
Probe 2; Willmer et al. 2005; Faber et al. 2005), VVDS (VIRMOS-VLT Deep Survey; Zucca et al. 2005) 

^We have adopted Bj — B = 0.15 and Svega = Sab- Uncertainties are as published and may not account for the contribution of 
large-scale structure. 

^We have not determined uncertainties for other surveys, but uncertainties for M{10~^'^h^ Mpc~^ mag~^) are probably comparable 
to those for M* 

*^Values of <x without uncertainties denote Schechter function fits where a was fixed. 
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Table 7 

S-BAND Luminosity Density of Red Galaxies 





All Li[mmo,-,i(K!,^ 


Luniinosit\" D< 

Mb < + 1.0 


-risifv /n(10^/i 

Mb <- 


Mb < -li - 1.0 


Mb : M" - i.'T 


a 


0.20 <z < 0.40 


8.08 ± 1.02 


7.05 ±0.89 


5.09 ±0.64 


1.70 ± 0.21 


0.50 ± 0.06 


-0.28 ±0.04 


0.40 <z < 0.60 


8.47 ± 0.75 


7.49 ± 0.67 


5.42 ± 0.48 


1.80 ± 0.16 


0.53 ± 0.05 


-0.28 ±0.05 


0.60 <z < 0.80 


9.18 ±0.56 


7.62 ± 0.47 


5.03 ±0.31 


1.44 ± 0.09 


0.39 ± 0.02 


-0.55 ± 0.06 


0.80 <z < 1.00 


10.7 ± 0.8 


9.25 ± 0.68 


6.37 ±0.47 


1.95 ± 0.14 


0.55 ± 0.04 


-0.43 ± 0.05 


0.20 < z < 0.40 


8.29 ± 1.05 


6.85 ± 0.86 


4.60 ± 0.58 


1.36 ± 0.17 


0.37 ± 0.05 


-0.50 


0.40 < z < 0.60 


8.57 ± 0.76 


7.18 ± 0.64 


4.82 ± 0.43 


1.42 ± 0.13 


0.39 ± 0.03 


-0.50 


0.60 < 2 < 0.80 


9.10 ± 0.56 


7.68 ± 0.47 


5.16 ± 0.32 


1.52 ± 0.09 


0.42 ± 0.03 


-0.50 


0.80 < z < 1.00 


10.8 ± 0.8 


9.19 ± 0.68 


6.17 ± 0.45 


1.82 ± 0.13 


0.50 ± 0.04 


-0.50 



^The Sun's absolute magnitude is Mb ^ 5.48 (Bessell, Castelli, & Plcz 1998). 

'^When integrating over portions of the Schechter function, the resulting js values can be very sensitive to the assumed or measured 
values of M* and a. For example, when a is fixed 3b{Mb < M* — 1)/ jsiMs < M* + 1) is a constant irrespective of the redshift or 

■^As discussed in §5.2, small errors in assumed or measured values of M* and a can produce large errors in the measured luminosity 
density of Mb < M* — 1.5 red galaxies. 
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